Apollonian Circle Packings: Geometry and Group Theory 

III. Higher Dimensions 



Ronald L. Graham 
Department of Computer Science and Engineering 
University of California at San Diego, La Jolla, CA 92093-0114 

Jeffrey C. Lagarias 
Department of Mathematics 
University of Michigan, Ann Arbor, MI 48109-1109 

Colin L. Mallows 
Avaya Labs, Basking Ridge, NJ 07920 

Allan R. Wilks 
AT&T Labs, Florham Park, NJ 07932-0971 

Catherine H. Yan ^ 
Department of Mathematics 
Texas A&M University, Cohege Station, TX 77843-3368 

(March 10, 2005 version) 

ABSTRACT 

This paper gives n-dimensional analogues of the Apollonian circle packings in parts I and II. 
Those papers considered circle packings described in terms of their Descartes configurations, 
which are sets of four mutually touching circles. They studied packings that had integrality 
properties in terms of the curvatures and centers of the circles. Here we consider collections of 
n-dimensional Descartes configurations, which consist of n -|- 2 mutually touching spheres. 

We work in the space M"^ of all n-dimensional oriented Descartes configurations parametrized 
in a coordinate system, ACC-coordinates, as those (n + 2) x (n -|- 2) real matrices W with 

W^Qd^„W = Clw,n where QD,n = x{ + ... + x^_^2 - ^(^^i H 1" Xu+t)"^ is the n-dimensional 

Descartes quadratic form, Qw,n = —8x1X2 + 2x3 -|- • • • -|- 2x^^2 5 ^"^^ QD,n and Qw,n are their 
corresponding symmetric matrices. On the parameter space Ai^ of ACC-matrices the group 
Aut{QD,n) acts on the left, and Aut{Qw,n) acts on the right. Both these groups are isomor- 
phic to the (n -|- 2)-dimensional Lorentz group 0(n + 1,1), and give two different "geometric" 
actions. The right action of Aut{Q]Vn) (essentially) corresponds to Mobius transformations 
acting on the underlying Euclidean space M" while the left action of Aut(Q£)^n) is defined only 
on the parameter space M"^- 

^Partially supported by NSF grants DMS-0070574, DMS-0245526 and a Sloan Fellowship. This author is 
also affiliated with Dalian University of Technology, China. 



We introduce n-dimensional analogues of the Apollonian group, the dual Apollonian group 
and the super- Apollonian group. These are finitely generated groups in Aut{QD,n), with the 
following integrality properties: the dual Apollonian group consists of integral matrices in all 
dimensions, while the other two consist of rational matrices, with denominators having prime 
divisors drawn from a finite set S depending on the dimension. We show that the Apollonian 
group and the dual Apollonian group are finitely presented, and are Coxeter groups. We define 
an Apollonian cluster ensemble to be any orbit under the Apollonian group, with similar notions 
for the other two groups. Wc determine in which dimensions there exist rational Apollonian 
cluster ensembles (all curvatures rational) and strongly rational Apollonian sphere ensembles 
(all ACC-coordinates rational). 

Keywords: Circle packings, Apollonian circles, Diophantine equations, Lorentz group, Coxeter 
group 

1. Introduction 

In part I we considered Apollonian circle packings, and observed that there exist such packings 
where the circles all had integer curvatures and the quantities (curvature) x (circle center) 
had integer entries. We gave an explanation for this phenomenon, in terms of the Descartes 
configurations in the packing. A Descartes configuration is a configuration of four mutually 
tangent circles, with all six tangency points distinct. We introduced the space A^d of all 
ordered, oriented Descartes configurations, parametrized in a coordinate system expressed in 
terms of the curvatures and centers of the circles in the packing. 

The explanation of the integrality properties was as follows. This configuration space of 
ordered, oriented Descartes configurations was shown to be a principal homogeneous space for 
the Lorentz group 0(3, 1), a six-dimensional real Lie group which we identified with the real 
automorphism group Aut{QD) of a quaternary quadratic form 

QDixi,X2,X3,X4) = xl+ xl+ xl + xl - ^{xi + X2 + X3 + ^4)^, 

which we called the Descartes form. We showed that the set of all Descartes configurations 
in the packing formed an orbit of a discrete subgroup A acting on this space, which was 
independent of the packing, and which we termed the Apollonian group. This group consisted 
of 4 X 4 integer matrices, and the integrality properties of the Apollonian group explained the 
occurrence of packings with integral curvatures and (curvature) x (circle center) data. If a single 
Descartes configuration in the packing has such integrality properties, then every Descartes 
configuration in the packing inherits such properties, and thus all the circles in the packing 
have such properties. 

Part I also introduced two other discrete subgroups of Aut^Qu), which had integrality 
properties and a geometric interpretation. These were the dual Apollonian group A'^ and the 
group generated by both the Apollonian group and the dual Apollonian group together, which 
we called the super-Apollonian group 

In this paper we study to what extent these results carry over to the n-dimensional case. We 
show that there exist n-dimensional analogues of the Apollonian group, dual Apollonian group 
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and super- Apollonian group. We study properties of their orbits on the parameter space 
of n-dimensional ordered, oriented Descartes configurations which we call Apollonian cluster 
ensembles. We determine to what extent properties of the two-dimensional case extend to the 
n-dimensional case. One immediate difference is that in dimensions n > 4 the orbits of the 
Apollonian group no longer correspond to sphere packings; the spheres in these configurations 
overlap. However the groups involved still have rationality properties, and we consider the 
question when there exist ensembles having all curvatures rational and/or having all data 
(curvature) X (sphere center) rational. 

In §2 we describe six properties that hold in the two-dimensional case, and then summarize 
what the main results of this paper say about generalizations of these to higher dimensions. 
Some of them generalize completely, others only in an infinite set of specific dimensions, and 
some are specific to dimension 2. 

In §3 we describe the parameter space A^jg, of all n-dimensional ordered, oriented Descartes 
configurations in a coordinate system, ACC-matrices, as those (n-|-2) x (n-|-2) real matrices W 
with W-^Q^^nW = Qw,n where QD,n = xl + ...+x'^_^_2 — ^{xi + ...+Xn+2)'^ is the (matrix of the) 
n-dimensional Descartes quadratic form and Qw,n = —8x1X2 + 2x3 + ... + 2x'^_^_2■ On the space 
Alp of ACC-matrices the group Aut{QD,n) acts on the left, and Aut(Qw,n) acts on the right. 
Both these groups are isomorphic to the (n -|- 2)-dimensional Lorentz group 0(n -|- 1, 1), and 
give two different "geometric" actions. The right action of Aut{Qw^n) (essentially) corresponds 
to Mobius transformations acting on the underlying Euclidean space while the left action 
of Aut(QD,n) is defined only on the parameter space. 

In §4 we define the n-dimensional analogues of the Apollonian group, the dual Apollonian 
group and the super-Apollonian group introduced in part I. As we just noted above, one 
immediate difference with the 2-dimensional case is that in dimensions n > 4 the orbits of the 
Apollonian group no longer correspond to sphere packings; the spheres in these configurations 
overlap. Furthermore, even viewed in the parameter space the action of the n-dimensional 
Apollonian group An is not discrete. However the action of the dual Apollonian group is 
discrete on Alp for all n > 2. To restore a discrete group action for the Apollonian group and 
the super-Apollonian group in certain higher dimensions, it suffices to view it as a diagonal 
action on M'^{M.) x flp j n-i -^b(Qp)) which AfJj(Qp) is a p-adic parameter space. This can 
be done in dimensions for which there arc rational points in the parameter space Atj^ = Alp(M). 

In §5 we determine presentations for the Apollonian group An in all dimensions n > 3. In 
dimension 3 there are extra relations of the form (SjSj)^ = I whenever i 7^ j. The group An 
is a hyperbolic Coxeter group in all dimensions n > 3. 

In §6 we investigate rationality properties of Apollonian cluster ensembles. We prove that 
strongly rational configurations exist if and only if the dimension n = 2/c^ or n = (2k + 1)^. 

In §7 we consider to what extent the duality operator D defined in dimension 2 has an 
n-dimensional analogue. In dimension 2 this operator takes a Descartes configuration P to a 
dual Descartes configuration D' consisting of the 4 circles each of which is orthogonal to all 
but one of the circles in the packing. We observe that this operation has a geometric analogue 
in dimension n > 3, which maps the configuration D to the set £^ of n + 2 spheres having the 
property of being orthogonal to all but one of the spheres in V. These spheres in 6 do not 
form a Descartes configuration, but have an equi-inclination property instead. In dimensions 
n > 3 this operation docs not have an algebraic interpretation in Aut(QD,n)- 

In §8 we make some concluding remarks, stating some unresolved questions. 
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In the Appendix we describe the action of the n-dimensional Mobius group Mob(n) on 
(ordered, oriented) Descartes configurations, encoded via an isomorphism, as the action of 

Aut{Qw,n)- 

The general framework of this paper was developed by the second author (JCL), who also 
did much of the writing. This paper is a revised and extended version of a preprint originally 
written in 2000. 

Acknowledgments. The authors are grateful for helpful comments from Andrew Odlyzko 
and N. J. A. Sloane. The authors thank the reviewer for incisive comments leading to reorga- 
nization of the paper. 

2. Main Results 

The object of this paper is to study to what extent the properties of Apollonian circle packings 
studied in parts I and II have n-dimensional analogues. The two-dimensional case had the 
following six features: 

(PI) [Parameter Space Property] The space Aio of all ordered, oriented Descartes con- 
figurations can be identified with the set of real intertwining matrices W under congruence 
between two rational quadratic forms in 4 variables, the Descartes form Qd and a certain form 
Qw, i-c. W-^Q/jW = Qw- Wc called the 4x4 matrix W the augmented curvature-center 
coordinates of the associated ordered, oriented Descartes configuration. The parameter space 
Mo is a real algebraic variety and a principal homogeneous space for the Lorentz group 0(3, 1). 

(P2) [Orbit Property] There exist three groups of matrices in Aut[Q£,) of determinants ±1, 
the Apollonian group A, the dual Apollonian group A'^ and the super- Apollonian group „4'^, 
which have geometrically characterizable actions on Descartes configurations. In particular, 
the set of Descartes configurations in an Apollonian circle packing comprise a single orbit of a 
single Descartes configuration under the action of the Apollonian group (disregarding orienta- 
tion and ordering). 

(P3) [Integer Matrix Property] The Apollonian group, dual Apollonian group and super- 
Apollonian group each consist of integer matrices. Thus all three groups are discrete subgroups 
of GL(4, R). They each have a discontinuous action on the parameter space M.^ of ordered, 
oriented Descartes configurations. 

(P4) [Coxeter Group Property] Each of the Apollonian group, dual Apollonian group, and 
super- Apollonian group have group presentations as hyperbolic Coxeter groups. 

(P5) [Integral Packing Property] There exist integer Apollonian circle packings, ones in 
which all Descartes configurations have integer curvatures. Furthermore there exist super- 
integral Descartes configurations, ones whose augmented curvature-center coordinate matrices 
W are integer matrices. There are super-integral Apollonian packings, ones in which ev- 
ery Descartes configuration has this property. The set of all super-integral ordered, oriented 
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Descartes configurations forms 672 orbits under tlie action of the super- Apollonian group. 

(P6) [Duality Operator] There exists a duality operator D in Aut{Qx>) given by a fixed 
matrix GL(4, M) with half-integer entries. It acts by conjugacy on the super- Apollonian group 
and gives an outer automorphism of order 2 of this group. This outer automorphism conjugates 
the Apollonian group to the dual Apollonian group. The duality operator has a geometric in- 
terpretation in terms of its action on individual Descartes configurations. 

The results of this paper generalize each of these features, which we consider in order. 
The parameter space property (PI) generalizes to all dimensions, as was shown by three of 
the authors in |2SI. Results in this direction were established earlier by Wilker [HSj, who 
used the term cluster for an ordered, but not oriented, Descartes configuration. There is a 
notion of augmented curvature-center coordinates (ACC-coordinates) for an ordered, oriented 
n-dimensional Descartes configuration, and a parameter space Alp of all such configurations, 
specified by a matrix condition that intertwines two quadratic forms in n -|- 2 variables under 
conjugacy, as was shown in [2S|- Furthermore the space Alp the structure of a principal 
homogeneous space for 0{n + 1, 1). We present these results in §3. 

The orbit property (P2) generalizes to all dimensions. There exists an n-dimensional version 
of the Apollonian group, the dual Apollonian group and the super- Apollonian group, with the 
same geometric interpretation of their actions on Descartes configurations. These are described 
in §4. We define an Apollonian cluster ensemble to be an orbit of a single (ordered, oriented) 
Descartes configuration under the action of the Apollonian group. However there is a sphere 
packing interpretation of the geometric packing corresponding to this orbit in dimension 3 
only, and not in dimensions n > 4. 

The integer matrix property (P3) partially generalizes to all dimensions, as we discuss in §4. 
The n-dimensional Apollonian group, dual Apollonian group and super- Apollonian group each 
consist of integral matrices in dimension 3, and they act discretely. The dual Apollonian group 
consists of integer matrices in all dimensions, and so is a discrete subgroup of GL{n + 2,M) and 
acts discontinuously on the parameter space Aip. However for n > 4 the Apollonian group and 
super-Apollonian group consist of rational matrices, with denominators whose prime factors 
all divide n — 1. The Apollonian group An is not a discrete subgroup of GL{n + 1,M) for 
n > 4 and does not have a discontinuous action on Al^ (Theorem 14. However we can 
restore discreteness by adding an action on some p-adic groups, for those primes p dividing 
n — 1. That is, the super-Apollonian group A^ can be embedded discretely by a diagonal 
action inside GL[n -|- 2, M) x | GL[n + 2, Qp) for all n. For certain dimensions n, those 
characterized in §6 as n = 2fc^ or n = {2k + 1)^ for /c > 1, one can get a discontinuous diagonal 
action of the super- Apollonian group on a parameter space Aip(M) x | „_i Aii[}(Qp)- For 
this it is necessary that the p-adic parameter spaces Aip(Qp) have enough points. Finally 
on the level of "packings" in the weak sense of spheres not crossing each other, this holds in 
dimensions 2 and 3 for all these groups, and holds for the dual Apollonian group A'^ for all 
n > 4 (Theorem 14. 3|) . but not for the other two groups. 

The hyperbolic Coxeter group property (P4) partially (and perhaps completely) generalizes 
to all dimensions. We show in §5 that for all n > 3 the n-dimensional Apollonian group is 
a hyperbolic Coxeter group (Theorems 15.11 and 15. 2|) . For n = 3 the group has an extra 
relation, which explains the existence of the structures "The Bowl of Integers" and "The 
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Hexlet" studied by Soddy |SJ, [SU and Gosset |T7j. It may be true that in all dimensions 
the dual Apollonian group and super-Apollonian group are also hyperbolic Coxeter groups, 
but we leave these as open questions. 

The integral packing properties in (P5) generalize to all dimensions in weakened forms. 
In §6 we observe that in all dimensions n > 3 there exist Apollonian cluster ensembles in 
which every sphere has a rational curvature, with denominators divisible only by a certain 
finite set of primes, those dividing n — 1 (Theorem 16. ip . We then consider the question in 
which dimensions n does there exist a super-rational Descartes configuration, i.e. one whose 
augmented curvature-center coordinate matrix Wx) is a rational matrix in A^jj. We prove this 
can be done if and only if the dimension n = or n = {2k + 1)^ for some k > 1 (Theorem 



The duality operation in (P6) giving a conjugacy between the Apollonian group and the 
dual Apollonian group, appears to exist algebraically only in dimension 2. In §7 we show 
that viewed as a geometric action, there does exist a natural "duality operator" acting on 
Descartes configurations in all dimensions n > 2. However in dimensions n > 3, the image of 
this operator is not a Descartes configuration, but instead is a collection of n + 2 spheres each 
of which intersects the others in a particular fixed angle 6n, which depends on the dimension 
n (Theorem rr^ . 

3. Descartes Configurations and Group Actions in M" 
3.1. Descartes Configurations and ACC-Coordinates 

We start with the generalization of the Descartes circle theorem to the n-dimensional case. An 
Descartes configuration in M" consists of (n+2) pairwise tangent (n— l)-spheres (5i, 52, • • • , Sn+2) 
in M"', with all points of tangency distinct. This result, which was termed the Soddy-Gossett 
theorem in after Soddy jHU] and Gossett U^l (although it was discovered earlier in the 
3-dimensional case), states that if the spheres have disjoint interiors then 



The Descartes circle theorem is the case n = 2. 

The Soddy-Gossett theorem holds for all Descartes configurations, including configurations 
where one sphere encloses the others, provided that we assign appropriate signs to the curva- 
tures, so that the configuration has a (total) orientation, as defined below. An oriented sphere 
is a sphere together with an assigned direction of unit normal vector, which can point inward 
or outward. If it has radius r then its oriented radius is r for an inward pointing normal and 
— r for an outward pointing normal. Its oriented curvature (or "signed curvature") is 1/r for 
an inward pointing normal and —1/r for an outward pointing normal. By convention, the 
interior of an oriented sphere is its interior for an inward pointing normal and its exterior for 
an outward pointing normal. An oriented Descartes configuration is a Descartes configuration 
in which the orientations of the spheres are compatible in the following sense: either (i) the 
interiors of all n + 2 oriented spheres are disjoint, or (ii) the interiors are disjoint when all ori- 
entations are reversed. Each Descartes configuration has exactly two compatible orientations 
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in this sense, one obtained from the other by reversing all orientations. The positive (total) 
orientation of a Descartes configuration is the one in which the sum of the signed curvatures 
is positive, while the negative (total) orientation is the one in which the sum of the curvatures 
is negative. One can check that the sum of the curvatures cannot be zero. 

Now let bi = ^ denote the (signed) curvature of the i-th sphere of an ordered, (totally) 
oriented Descartes configuration , and let b = (61, 6n+2)- The geometry of such a Descartes 
configurations is encoded in the signed curvature vector b. In the positively oriented case, 
where Yl^=i > 0' the following holds: (i) all of 61, 62, • • • , &n+2 are positive; (ii) n + 1 

are positive and one is negative; (iii) n + 1 are positive and one is zero; or (iv) n are positive 
and equal and the other two are zero. These four cases correspond respectively to the following 
configurations of mutually tangent spheres: (i) n + 1 spheres, with another in the curvilinear 
simplex that they enclose; (ii) n + 1 spheres inscribed inside another larger sphere; (iii) n + 1 
spheres with one hyperplane (the (n + 2)-nd "sphere"), tangent to each of them; (iv) n equal 
spheres with two common parallel tangent planes. 

We can reformulate the Soddy-Gossett theorem in matrix terms. Let b-i = ^ denote the 
(signed) curvature of the i-th sphere of an ordered, oriented Descartes configuration , and let 
b = (61, 6„+2), then p. If) becomes 

b^QD,„b = 0, (3.2) 

in which Q_D,n is the symmetric matrix associated to the n-dimensional Descartes quadratic 
form, defined below. 



Definition 3.1. The n-dimensional Descartes quadratic form QD,n is the quadratic form in 
n + 2 variables whose associated symmetric matrix is 



QD,n '■— In+2 



1 



n 



Ln+2-'-n+2) 



in which 1 = (1, 1, 1)"^ is a column vector of length n + 2. 
Here the original Descartes quadratic form is 



Q 



D,2 



(3.3) 



(3.4) 



In PSl three of the authors of this paper showed that there exists a parametrization of 
the set of all ordered, oriented Descartes configurations in M", using a coordinate system 
involving the curvatures and centers of the spheres, which appears as Theorem 13 . 1 1 below. 



Definition 3.2. Given an oriented sphere S in M", its curvature- center coordinates consist 
of the (n + l)-vector 

m(S') = (6,6x1,... , ten), (3.5) 
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in which b is the signed curvature of S (assumed nonzero) and x(iS') = x = {xi,X2, ■ ■ ■ , Xn) is its 
center. For the degenerate case of an oriented hyperplane H, its curvature-center coordinates 
m{H) are defined to be 

m(5) = (0,h), (3.6) 

where h := (/ii,/i2, . . . ,hn) is the unit normal vector that gives the orientation of the hyper- 
plane. 

To see the origin of this definition in the degenerate case, let the point of H closest to the 
origin be z = Ah for some real value A. For t > |A|, let St be the oriented sphere of radius 
t centered at (t + A)h, which has center in direction h from the origin and contains z. As 
t OC' the oriented spheres St clearly converge geometrically to the oriented hyperplane H, 
and m(5t) = (i, (1 + f )h) -> ni{H) = (0, h). 

Curvature-center coordinates are not quite a global coordinate system, because they do not 
always uniquely specify an oriented sphere. Given m G M"+^, if its first coordinate b is nonzero 
then there exists a unique sphere having m = m(iS'). But if 6 = 0, the hyperplane case, there 
is a hyperplane if and only if ^ /if = 1, and in that case there is a pencil of hyperplanes that 
have the given value m, which differ from each other by a translation. 

We obtain a global coordinate system for spheres by adding an additional coordinate. This 
coordinate incorporates information about the sphere S obtained from S by inversion in the 
unit sphere. In n-dimensional Euclidean space, the operation of inversion in the unit sphere 
replaces the point x by x/|x|^, where |xp = Yl^=i^^- Consider a general oriented sphere S 
with center x and oriented radius r. Then inversion in the unit sphere takes S to the sphere 
S with center x = x/(|xp — r^) and signed radius f = r/(|xp — r^). If |x|^ > r^, then S has 
the same sign as S. In all cases, 

X X 

r r 

and 

= r 

r 

Definition 3.3. Given an oriented sphere S in M", its augmented curvature-center coordinates 
(or ACC- coordinates) of S are given by the (n -|- 2)-vector 

w(S) := (5, 6, 6x1, ... , bxn) = {b, m), (3.9) 

in which b = b{S) is the curvature of the sphere or hyperplane S obtained by inversion of S in 
the unit sphere, and the entries of m are its curvature-center coordinates. For hyperplanes we 
define 

wiH) := ib,0,hi,..., hn) = (b, m), (3.10) 

where b is the oriented curvature of the sphere or hyperplane H obtained by inversion of H in 
the unit sphere. 

7 



(3.7) 



(3.8) 



Augmented curvature-center coordinates provide a global coordinate system for oriented 
spheres: no two distinct oriented spheres have the same coordinates. The only case to resolve 
is when 5 is a hyperplane, i.e., 6 = 0. The relation (|3.7|) shows that {b, bxi, . . . , hxn) are the 
curvature-center coordinates of S, and if 6 7^ 0, this uniquely determines S; inversion in the unit 
circle then determines S. In the remaining case, 6 = 6 = and S = S is the unique hyperplane 
passing through the origin whose unit normal is given by the remaining coordinates. 

Given a collection {Si, S2, ■ ■ ■ , Sn+2) of n + 2 oriented spheres (possibly hyperplanes) in M", 
the augmented matrix W associated with it is the (n + 2) x (n + 2) matrix whose j-th row has 
entries given by the augmented curvature-center coordinates w(5j) of the j-th sphere. 

To state the next result, we introduce another quadratic form. 

Definition 3.4. The (n- dimensional) Wilker quadratic form Qw,n is the (n + 2)-variable 
quadratic form given by the symmetric (n + 2) x (n + 2) matrix 



Qw,n ■~- 



0-4 
-4 
2I„ 



(3.11) 



This name is made in honor of J. B. Wilker [HS], who introduced in spherical geometry a 
coordinate system analogous to augmented curvature-center coordinates, see [23 §2 p. 388-390 
and §9]. However he did not formulate any result explicitly exhibiting a quadratic form like 
Q,w,n'i see the remark on [2S1 P- 349 ]. 

Theorem 3.1. [Augmented Euclidean Descartes Theorem] The augmented matrix W = Wx) 
of an oriented Descartes configuration D of n + 2 spheres {Si : 1 < i < n + 2} in satisfies 






-4 





-4 














21 



(3.12) 



Conversely, any real solution W to 
oriented Descartes configuration D. 



is the augmented matrix Wxi of a unique ordered, 



Proof. This is proved as Theorem 3.3 in |^. ■ 



Theorem 13 . 1 1 st at es that the augmented curvature-center coordinates of an ordered, oriented 
Descartes configuration give an intertwining map between the Descartes form and the Wilker 
form. We note that the Soddy-Gossett theorem is a special case of Theorem 13.11 encoded as 
the (2, 2) entry of the matrix W^Q/)^„W. 

Both the Descartes quadratic form and the Wilker quadratic form are equivalent over the 
real numbers to the Lorentzian quadratic form 

QcA^) ■= -^0 + h x^+i, (3.13) 

see §3.2. This quadratic form has a large group of real automorphisms under congruence 
Aut{Qc,n) = {U G GL(n,M) : U^Q£,nU = Qc,n}, 
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which is the Lorentz group 0(n + l, 1). In consequence both the Descartes quadratic form QD,n 
and Wilker quadratic form Qw,n have automorphism groups under (real) congruence which 
are conjugates of the Lorentz group. 

The Descartes form is not only equivalent to the Lorentz form over the real numbers, but 
sometimes over the rational numbers. In dimension 2, the Descartes and Lorentz forms are 
rationally equivalent, where one has 



Jo"Q-D,2Jo, 



(3.14) 



with 



1 1 

-1 -1 

1 -1 

-1 1 



In dimension 2 the Wilker form Qw,2 is also rationally equivalent to both the Descartes form 
and Lorentz form, as shown in part I, §3.1. 

In higher dimensions these three forms are not always rationally equivalent. A necessary 
condition for rational equivalence of two quadratic forms is that their determinants differ by 
a rational square. We have det(Q£) = —1, det(Q£)) = — ^ and det(Qiy) = —2""'"^. It follows 
that a necessary condition for rational equivalence of the Descartes and Lorentz forms in 
dimension n is that n = 2fe^; for the Wilker and Lorentz forms in dimension n that n = 2k; 
and for the Descartes and Wilker forms in dimension n that n = k'^ for odd k or n = 2k'^ for 
even k. All three of these necessary conditions hold if and only if n = 2k'^. 

In §6 we show that the last condition is sufficient for equivalence of the Descartes and 
Wilker forms. 



3.2. Mobius and Lorentz group actions 

The augmented Euclidean Descartes theorem yields two group actions on the space of Descartes 
configurations. The group Aut{QD^n) acts on the left and the group Aut{Qw,n) acts on the 
right, as 

Wt)^UWt)V^, with V e Aut{QD,n), ^ Aut{Qw,n)- 

The two group actions obviously commute with each other. The following result generalizes 
the two-dimensional case, where the Lorentzian form Qc,2 is given in (|3.14|) . 

Theorem 3.2. (1) The groups Aut{QD,n) and Aut{Qw,n) ctre conjugate over the real numbers 
to Aut{Qc,n) = 0{n + 1,1). 

(2) The group Aut{Qo,n) acts transitively on the left on the parameter space of all 
ordered, oriented Descartes configurations. Given two such Descartes configurations T> and T> 
there exists a unique U G Aut{Q£),n) such that UW^i = W^)'. 

(3) The group Aut{Qw) acts transitively on the right on the space of all ordered, oriented 
Descartes configurations A^jj- Given two such Descartes configurations T> and T> there exists 
a unique V G Aut{Qw,n) such that W©V"i = Wp/. 
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Remark. This result allows one to define both a left and right 0{n + 1, 1) action on the 
space M.^, which depends on the choice of conjugacy made in (1). Then (2) and (3) show 
both these actions are transitive and have trivial stabilizer. This gives the structure of a 
principal homogeoneous space (or torsor) for 0(n + 1, 1), for each of these actions. 



Proof. Part (1) follows for Qw,n on taking Qc,n = Z-^Qvi/,nZ with 



1 

71 



r 1 _i 

? ? 



2 2 

1 







It then follows for QD,n because it is conjugate to Qw,n by Theorem 13.11 using for any 
fixed Descartes configuration. 

Parts (2) and (3) follow immediately from (1), because Aut{Qc,n) = 0{n + 1, 1) acts tran- 
sitively. ■ 



Since the spheres in a Descartes configuration appear as the rows in the augmented matrix 
Wx) of an oriented Descartes configuration, the action on the right by elements of Aut{Q\Y^ri) 
maps spheres to spheres. This action can essentially be identified with the Mobius group of 
linear fractional transformations acting on the one-point compactification M" of M". More 
precisely, it corresponds to a direct product of the Mobius group with {I, —I}, because the 
Mobius group preserves total orientation of Descartes configurations. A precise description of 
the isomorphism is given in the Appendix. (The case n = 2 was treated in Appendix A of part 

I-) 

The action on the left, by Aut{QD,n), mixes together the different spheres in the original 
Descartes configuration, and does not make sense as an action on individual spheres. This 
group action is intrinsically associated to the "("+^) -dimensional (real) parameter space 
of oriented Descartes configurations. 



4. Apollonian Groups and Apollonian Cluster Ensembles 

In parts I and II we studied Apollonian circle packings in terms of the Descartes configura- 
tions they contain. We showed they consisted of a single orbit of a discrete subgroup of the 
automorphism group Aiit((5/3 2) of the Descartes quadratic form Qd,2- We also introduced 
there another discrete subgroup of Aut{QD^2), the dual Apollonian group, and in addition, 
the super-Apollonian group, which is the subgroup of Aut((5D,2) generated by the Apollonian 
group and the dual Apollonian group together. 

We show there are analogues of these three groups in all dimensions n > 3, which are 
subgroups of Aut^Qu^n) consisting of rational matrices. We call these the n-dimensional 
Apollonian group, dual Apollonian group, and super-Apollonian group. We then define an 
Apollonian cluster ensemble to be an orbit of the Apollonian group. This provides a gener- 
alization of Apollonian packing to all dimensions, though it turns out not to correspond to a 
sphere-packing in dimensions n > 4. 
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4.1. n-Dimensional Apollonian Group 

The n-dimensional Apollonian group An = (Si, S2, S„+2) consists of (n + 2) x (n + 2) 
matrices with 

S, := I + -^ejl^ - -^ejej, (4.1) 
n — 1 n — \ ■' 

where e-^ is the j-th unit coordinate (column) vector and 1 = ei + • • • + e„+2 = (Ijl;---)!)"^- 
That is, Sj is the identity matrix in all rows but the j-th row, and there has —1 on the diagonal, 
and all off diagonal elements equal to It is straightforward to check that An C PLU.t{QD,n)- 
The relations = I are evident. 

The algebraic action of the operator Sj on the augmented curvature center coordinates Wx> 
of a Descartes configuration P is to take it to SjWx>, which is Wj?/ for some V . The geometric 
interpretation of this action is to fix all (oriented) spheres in the Descartes configuration T> 
except the j-th sphere, and to replace that sphere with the unique other sphere that is tangent 
to the remaining n + 1 spheres, assigned an appropriate orientation; this is T>' . For a fixed 
Descartes configuration P this operation can also be realized by a Mobius transformation that 
is an inversion with respect to the sphere that passes through the "^-"^^"^^ tangency points of 
the remaining n + \ spheres in the Descartes configuration. (The existence of such a sphere is 
demonstrated in Theorem 17. ll ) 

It is apparent that An is a group of integer matrices for n = 2 and 3, while for n > 4 it 
is not always integral, consisting of rational matrices whose denominators contain only powers 
of primes that divide n — 1. 

Theorem 4.1. The Apollonian group An is a discrete subgroup of GL(n + 2,]R) for n = 2 
and n = 3. It is not a discrete subgroup of GL{n + 2,M) for all n > 4. 

Proof. The group An is discrete in dimensions 2 and 3 since it is then a subgroup of 
GL{n + l,Z). 

We consider the element S1S2, and show for n > 4 that the set {(8182)'^ : A; > 1 does not 
contain the identity matrix, but the closure of this set does contain the identity matrix. The 
(n + 1) X (n + 1) matrix 8182 is a product of two reflections, so it has determinant 1 and has 
n of its eigenvalues equal to 1. Its first two rows are (a^ — 1, — a„, + an,afi + a„, ...) and 
{on, — 1, a„, a„, ...) in which a„ = ^^r:]-' ™^ other rows are those of the identity matrix. Its 
two non-unit eigenvalues are e^*^" with cos ^ = so 0„ = 2 cos^"'^(^^). Since these two 
eigenvalues are distinct, it follows that 8182 is diagonalizable. Now := {(8182)'^ : A; > 1} 
contains diagonalizable elements with all eigenvalues arbitrarily close to 1, so it contains the 
identity matrix as a limit point. When n > 4 it is well-known that On is an irrational multiple 
of vr, whence S„ does not contain the identity matrix, which shows that the Apollonian group 
A n is not a discrete subgroup of GL{n + 1,M). ■ 

Because the Apollonian group is arithmetic one can view it as a discrete group if one 
considers its actions on certain p-adic groups GL{n + 2,Qp) for p\n — 1. Let n* = n — 1 if 
n is even, and if n is odd. Then one can establish that the Apollonian group An is a 
discrete group when embedded diagonally in GL{n -|- 2, M) x Y\p | „. GL{n + 2, Qp). Using this 
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structure, one may obtain a discrete (i.e. discontinuous) diagonal action of the Apollonian 
group on the space M.^(M.) x Y[p | „* ■M.jiiQp), whenever all of the parameter spaces A^p(Qp) 
are nonempty. Here the spaces M.'^{Qp) are defined as the p-adic solutions to the conditions 
in (|3.12|) . The criterion in §6.2, shows that A^p(M) contains rational solutions, hence p-adic 
solutions for all p, when n = 2k'^ or (2k + 1)^, for some fc > 1. 

4.2. n-Dimensional Dual Apollonian Group 

The n- dimensional dual Apollonian group A:^ = {Sj^, S^, S;^_,_2) consists of (n + 2) x (n + 2) 
matrices with 

Sf := I + 21eJ - AejeJ. (4.2) 

That is, Sj~ is the identity matrix in all columns except the j-th column, where it has —1 as a 
diagonal entry and 2 for all off-diagonal entries. The group is a group of integer matrices 
in all dimensions n. It is straightforward to check that A^ C Aut((5£)^„). It is evident that 
(Sj")'^ = / holds for all j. The geometric interpretation of the operation Sj" on a Descartes 
configuration D is that it encodes inversion with respect to the i-th sphere of that configuration. 

Theorem 4.2. The dual Apollonian group A:^ is a discrete subgroup of GL{n + 2,R) for all 
n > 2. 

Proof. This holds since A^ is a subgroup of the discrete group GL{n + 2,Z). ■ 

Orbits of this group, acting on Descartes configurations, retain a "packing" property at the 
level of individual spheres. 

Theorem 4.3. For all n > 2, an orbit of the dual Apollonian group acting on a single 
Descartes configuration gives a "packing" of spheres in the weak sense that no two spheres in 
distinct Descartes configurations of the orbit cross each other, i.e. any two such spheres either 
coincide, or are disjoint, or are tangent. 

Proof. This can be seen geometrically by constructing the dual packing starting from a 
single Descartes configuration D. The spheres in the resulting "packing" are all nested inside 
spheres of P, and we call the level of a sphere its depth of nesting inside some sphere of P, 
the spheres in V being assigned level 0. One proceeds in stages, where at stage k — 1 one 
has the set of Descartes configurations {S^_ S^_ • • • S^Wx)}, in which each ij ^ ij-i- We 
assert as an induction hypothesis, that each stage k — 1 Descartes configuration has all but 
one of its spheres at level k — 1, with one sphere at level k — 2. This holds for the base case 
/c = 2 by inspection. At stage k each multiplication by a generator takes a particular Descartes 
configuration at stage k — 1 and maps (n + 1) of its spheres inside one level k — 1 sphere of the 
configuration. (Here the condition ik ^ ik-i is used.) The new Descartes configuration then 
consists of (n + 1) spheres nested to depth k inside V, contained in one outer sphere nested to 
depth k — 1. Furthermore each level k Descartes configuration is nested inside a unique level 
k — 1-sphere, of which there are (n + 2) • (n + 1)^~^ choices. It follows that the new level k 
spheres cannot cross any spheres at any levels up to k — 1, and they also cannot cross any 
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other level k spheres because they are either inside different level k — 1 spheres, or if they are 
in the same level k — 1 sphere, then they form part of a single Descartes configuration. Thus 
the non-crossing property holds at depth k, and the result follows by induction on A;. ■ 

4.3. n-Dimensional Super-Apollonian Group 

The n- dimensional super- Apollonian group is the group generated by An and A^, so that 

This group consists of integer matrices when n = 2 or 3, and of rational matrices otherwise. In 
particular A^ is a discrete subgroup of GL{n + 2,R) for n = 2 or n = 3, and is not a discrete 
subgroup for n > 4. The property that some spheres in different Descartes configurations 
intersect non-tangentially when n > 4 is inherited from the action of the Apollonian group. 

We do not know if the super- Apollonian group is a Coxeter group. We do know that its 
generators satisfy the Coxeter relations given in the following lemma. It seems plausible that 
for n > 4 these are a generating set of relations; if so, then Af^ would be a hyperbolic Coxeter 
group. 

Theorem 4.4. The super- Apollonian group Af^ = (Si, S2, • • • S„+2j S]*-, S2", S;[-_|_2) has gen- 
erators satisfying the relations. 

S] = I and {Sf f = /, 1 < j < n + 2. (4.3) 

and 

SjSi = S-tSj when j ^ k. (4.4) 

Proof. This is a direct calculation from the definition. To appreciate 1)4. 4() one can consider 
more generally (n + 2) x (n -|- 2) matrices Sj(A) whose entries are the identity matrix except 
in the j-th row, where they are —1 on the diagonal and A off the diagonal. Then Sj(A)^ = / 
holds for 1 < j < n, but the extra relations 1)4.4(1 hold only when A = (Compare the 

(j, /c)-th entry of the products SjS^ and S^S'j). Note however that Sj(A) S Aut{QD,n) if and 
only if A = ;^. . 

As mentioned above for the Apollonian group, there is a discrete action of the super- 
Apollonian group in certain dimensions, given by its diagonal action of A^ on the product 
parameter space -A^d(IK) xflp | n* -^d(Qp)) provided that all the spaces M^{Qp) are nonempty. 
As shown in §6.2, this will be the case when n = 2k'^ or {2k + 1)^, for some k > 1. 

4.4. Apollonian Cluster Ensembles 

In the two-dimensional case Apollonian packings could be described as (a) collections of circles, 
or (b) the orbit of four circles under a certain discrete group of Mobius transformations or (c) 
the Descartes configurations given by an orbit of the Apollonian group. The first two of these 
notions do not generalize to all dimensions, but version (c) does. 
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Definition 4.1. An Apollonian cluster ensemble in n dimensions is defined to be the orbit 
of the Apollonian group An of a given ordered, oriented Descartes configuration Dq in the 
parameter space M.^. 

This object can be viewed as a discrete set in the parameter space Ai'^ of all (ordered, 
oriented) Descartes configurations in dimension n. In dimension n = 3 the spheres in an Apol- 
lonian cluster ensemble yield a sphere packing. However in dimensions n > 4 the spheres in 
any such ensemble overlap and no longer correspond to a packing, as shown by Lemma [4.1l Be- 
cause the Apollonian group consists of rational matrices, we can ask if there are n-dimensional 
Apollonian cluster ensembles with rationality properties, either of their curvatures or of their 
full augmented curvature-center coordinate matrices Wx). 

Similarly one can define a dual Apollonian cluster ensemble to be an orbit of the dual 
Apollonian group of a single Descartes configuration. It is necessarily a discrete set, because 
An has an integral structure, so the identity matrix is isolated in the group A^- We can ask 
about rationality or integrality properties of orbits of this group. 

A super- Apollonian cluster ensemble is an orbit of the super- Apollonian group A^. We can 
ask about rationality properties of orbits of this group, and integrality properties in dimensions 
n = 2 and n = 3. For n = 2 these were answered in part II, and for n = 3 we show in §6 that 
integral do not exist. For n > 4 these orbits are not discrete, but in certain dimensions 
can have a rational structure inherited from the super-Apollonian group, see §6. 

5. Presentation for n-Dimensional Apollonian Group 

In we obtained a presentation for the super-Apollonian group dimension 2, which 

established that it was a hyperbolic Coxeter group. Theor em l4 . 41 gave some non-trivial relations 
in the super-Apollonian group A^, all of Coxeter type. There are additional relations, at least 
in dimension 3, see Theorem 15. II It may well be that A^ is a hyperbolic Coxeter group for all 
n > 3. However it seems a complicated problem to determine a presentation of A^ in general, 
and here we establish a more limited result. 

We determine a presentation for the Apollonian group An in dimensions n > 3, given in 
the next two results. In all cases it is a hyperbolic Coxeter group. 

Theorem 5.1. For dimension n = 3 the Apollonian group A3 = (81,82,83,84,85) has the 
relations = I for 1 < j < 5 and the additional relations 



(SjSk)^ = I, when j ^ k. 



(5.1) 



These are a generating set of relations, and A3 is a hyperbolic Coxeter group. 



Proof. Recall that 



81 



-11111 
10 
10 
1 
1 



(5.2) 
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and the other Sj are permutations of the j-th and first rows and columns of Si. It is easy 
to check that the generators satisfy ah the relations given above; in what follows we denote 
this set of relations TZ. We note that relations H5.1() can be written in either of the forms 
(SjSfc)^ = SfcSj or SjSfcSj = SfcSjSfc. 

We now consider a general word U = U1U2 • • • U„ with each JJj = Si-, where inverses of 
generators are eliminated by relations S^^ = Sj. A subword of U is any word UjUj+i • • • U^. 
We call a word reduced if it has the following properties. 

(i) 2-reduced. It contains no subword of form SjSj. 

(ii) B-reduced. It contains no subword of form V1V2 • • - ^21x1 in which Vi = V3, = 
V2j+3 for 1 < J < m - 2, and V2m-2 = ^2m- 

Conditions (i) and (ii) together allow the refinement of (ii) to assert Vi = V3 7^ V2, 
V2j = V2j+3 / V2j+2, V2m-2 = ^2m / ^2m-^i- Also the case m = 2 rules out words 
SjSfcSjSfc. The definition also implies: every subword of a i?-reduced word is also i?-reduced. 

We assert that every non-reduced word can be simplified to a reduced word of a shorter 
length, using the relations IZ. If (i) is violated, then using the relation = / we may 
replace the subword by the empty word, shortening it by 2. If (ii) holds then we can replace 
the subword by V2V3 • • • V2m-i, decreasing its length by two. This reduction is done by a 
sequence of replacements of of the form SjSfcSj by SfcSjSfe applied successively at positions 
1,3, ...,2m — 3, which ends with V2 at the beginning and ^ 2m-i^ 21^1^ 2m at the end. Then 
the final letters V2mV2rri are deleted by rule (i), achieving a shortening by 2. Starting with 
any word U and applying this reduction process, we eventually arrive at a reduced word or 
the empty word. 

The theorem is equivalent to the assertion that each nonempty reduced word is not the 
identity. For if the relations IZ did not generate all relations, there would exist another relation, 
necessarily forcing some nonempty reduced word to be the identity. 

We introduce some (matrix) invariants associated to a word in the generators. Let 
be the j-ih. unit (column) vector, i.e. the j-th column of the 5x5 identity matrix, and set 
I5 = (1,1,1,1,1)^. We define 

C7,(U) :=eJUl5, (5.3) 

the sum of entries in row j of U, and 

S(U) := ll^Ul5, (5.4) 

the sum of all entries of U, which we call its mass. 

We compute the effect on these invariants of multiplying by a generator. The matrix SjU 
has all rows but the j-th row the same as U, with its j-th row equal to the sum of all rows of 
U minus twice its j-th row, which yields 

(7j{Sj\J) = S(U) - 2f7j(U), (5.5) 

and cTfc(SjU) = o"fc(U) if k ^ j. It also gives 

S(SjU) = 2S(U) - 3aj(U). (5.6) 

Now, we define 

6j{U) := S(SjU) - S(U). (5.7) 
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This quantity measures the increase in the total mass of a matrix when we multiply by Sj. 
Two key properties of this measure are that if j ^ k then 

(5,-(SfeU) = 5,-(U) + 5fc(U), (5.8) 

while if j = A: then 

6j{SjV) = -6j(U). (5.9) 

To verify these, first observe that 

6j{U) = (2S(U) - 3aj(U)) - E(U) = E(U) - 3aj(U). 



Therefore 



5,(SfcU) = S(SfcU)-3a,(S,,U) 

= 2S(U)-3(Tfc(U)-3CTj(U) 

= S (U) - 3aj- (U) + S (U) - 3afc (U) 

= Sj{V) + 6k{\J), 

which gives ()5.8|) . Similarly, using l\5.7\i . 

SjiSjV) = S(SjSjU) - S(SjU) 
= -(S(S,U) - S(U)) 

= -5.(u). 

One consequence of these two properties is that 

5j{SkU) = SkiSjV) if j^k. (5.10) 

Another consequence is 

5,(SfeS,U) = 6k{SjU) + S.iSjV) = (5,(U) + 6kiU)) - 6j{lJ) = 6k{lJ). (5.11) 

We assert that all nonempty reduced words U = SjU' have 

6j(U') > 0, (5.12) 

so that I1(U) > S(U'). If this is proved, then since U' is also reduced, we obtain by induction 
on the length of U that 

S(U) > S(U„) = 7 > 5](/) = 5, 

so that U cannot be the identity matrix, and the theorem follows. 

We establish (|5.12|) by induction on the length n of U = U1U2 • • • U„ = UiU', where we 
suppose Ui = Sj. In the base case n = 1, we consider that U' = I, and we then have 

S(Sj) = 7 > S(/) = 5, 
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completing the base case. Now suppose n > 2 and that the induction hypothesis holds up to 
n — 1. We write U = SjSfcU" with U" of length n — 2, noting that j ^ k. We are to show 
Sj{U') > 0. Now (EIHl) gives 

6j{V') = SjiSkV") = 6,{V") + 4(U"). (5.13) 

The induction hypothesis gives 6k(U") > 0. If SjU" is reduced, then (5j(U") > by the 
induction hypothesis and we are done. So suppose SjJJ" is not reduced. Since U" is reduced, 
a non-reduced subword in it must be an initial segment, which either fails to be 2-reduced or 
-B-reduced. 

Suppose first that S^U" is not i?-reduced, having an initial i?-word Vi • • • V2m which 
begins SjSiSj so that SjU" = SjSiSjY, for some I ^ j. We also have / / k since U' = SfcU" 
is 2-reduced. Now we have, using (|5.1()|) . 

6jiV') = 6j{SkV") = 6k{SjV"). (5.14) 

Applying the i?-reduction procedure to the initial segment Vi • • • \^2m simplifies the word SjU" 
to a word S^SjV that is shorter by two letters than SjJJ" , but is equal to it as a matrix. It 
follows using (|5.14() that 

5k{SiSjV') = SkiSjV") = 6j{V'). (5.15) 

We now assert that S^Sj-V itself is reduced. To see this, write S^Sj-V = V2V3 • • • V2m~i V' 
with V" = SfStV". Since U" = S^V was reduced, any further reduction involves a sub- 
word X that contains V2m-2Sr. Now 7^ V2m or else U" includes V2mSr = Y2m^2m 
contradicting U" being reduced. Next 7^ Y2m~i, for otherwise U" would contain the 
subword y^2m-2^2m-i^2my^2m~i which is a i3-word since \^2m = ^2m-2, contradicting U" 
being reduced. It follows that \'2m-2^m-iSr = ^2m^2m-iSr consists of three distinct letters 
(generators). Thus any non-reduced subword X = X'V2mV2m~iSrX" inside S^SjV must be 
a i?-word. But then U" contains the i?-word \'2rn^2m-i^2mSr^" , contradicting U" being 
reduced. We conclude that S;SjV' is reduced. We deduce that SfcS^SjV is reduced, because 
S/SjV is reduced and differs from both of its first two letters. The induction hypothesis 
now applies to give 5fc(S;SjV') > 0, which with ()5.15|) gives Sj{\J') > 0. 

The remaining case is that when SjU" is not 2-reduced. Then U" = SjV, and U = 
SjSfcSjV. We have then, by that 

If SfcV is reduced, then the induction hypothesis gives (5fe(V) > 0, which gives the desired 
result. If it is not reduced, then it is either not 2-reduced or S-reduced. If it is not 2-reduced, 
then we have V = S^V, in which case the original word U = SjSfcSjSfcV is not i?-reduced, 
contradicting the hypothesis. If S^V is not S-reduced, then we have V = S/S^V with I 7^ k, 
and U = SjSfcSjS/SfcV which shows that I ^ j- Now we have by ()5.11|) again that 

If SfcV is reduced then the induction hypothesis gives Si(V') > which implies that Sj{\J') > 
as desired. If it is not 2-reduced then V = S^V" and this contradicts U being reduced. There 
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remains the case where V' is not i?-reduced. Then we get V = SmS^V", with 5a;(V') = 
(5m, (V"). We can continue recursively in this way until the entire word U is used up and there 
are not enough letters to have any i?-reduced word. We then obtain 

6j{\J') = 6k{Y) = 6i(V) = --- = 5„(VW) > 0. 

This completes the induction step, proving H5.12() . and the theorem follows. ■ 

Theorem 5.2. For dimension n > 4 the Apollonian group An = (Si, S2, S„-|-2) 'is a hy- 
perbolic Coxeter group whose only relations are 

S]=ln+2 for l<j<n + 2. 



Proof. Suppose n > 4. We must show that no nonempty product of the Sj's, with distinct 
adjacent terms, is equal to the identity matrix In+2- 

Let Qn = and let Aj = — In+2) so that Sj = In+2 + 2Aj. Then Aj has all entries 
zero except for the jth row, in which all entries are gn except for the jth element, which is — 1. 
Let 

A = gnln+2^n+2 ~ n^n+2 

be the matrix with —1 on the main diagonal and gn elsewhere, so that its {i,j) element is 

/■ -N r if i = j, 

9n, if ^ / J- 

Then Aj and A both have jth row (a(j, 1), . . . , a(j, n + 2)). We note that A is nonsingular. 

Suppose that k >2, that 1 < j'l, • • • , j'fc < n + 2, that no two consecutive jV's are equal, and 
that S is the product 5 = Sj-^Sj^ • • • Sj^. We assume that S = In+2 and derive a contradiction. 
We have: 

S - In+2 = Sj^ Sj2 • • • Sj^^ - In+2 

= (In+2 + 2 Aj, ) • • • (I„,+2 + 2 Aj, ) - In+2 

= 2 5; A,, +4 A,,A,,+8 A,vA,,A,, + ... + 2%,...A,,. 

l<r<fc l<r<s<fc l<r<s<t<k 

The Aj's multiply more simply than the Sj's, compensating for the more complicated 
expression involving them. A product like SAj^Aj^Aj^, which we call an A-product, has one 
nonzero row, row jr, which is equal to row jt of A multiplied by 8a{jr, js)0'{js, jt)- The terms 
of the latter product are —1 01 g — n, depending on whether or not jr = js and js = jt- We 
call the scalar 8a{jr, js)0'{js, jt) the a-product corresponding to the A-product SAj^Aj^Aj^. In 
general, an a-product has the form ±2'^gn,a > /3. 

The jth. row of S — In+2 gets a contribution from each A-product that starts with Aj. It 
is a linear combination of rows of A determined by the last subscript in each contributing A- 
product. Because this linear combination is zero, each of its coefficients must be zero, since A is 
full rank. Therefore, for each t and u in {ji, ■■■,jk}, the sum of all the A-products corresponding 
to A-products that start with Af and end with A^ must be zero. This sum is a polynomial 
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in the variable g = gn, call it Pt^u{g)- Because A-products have the form :iz2"gf^,a > P, the 
coefficient of 5^ in this polynomial is divisible by 2*+^, so we may write 

fe-i 

i=0 

for integers Ct^u,i- 

Now consider Pji,j,.ig)- This has degree k — I. In fact, the only A-product that can 
contribute an a-product of degree fc — 1 is • • • A^^, and since successive A^^'s are distinct, 

the corresponding a-product is, 2'^g^~^, i.e., Cj^j^^^k-i = 1- Therefore, by the Rational Root 
Theorem, since gn = is a root of Pji,j^.{g), we must either have n — 1 dividing 2^^, or 
n = 2"* + 1 for some m > 1, in which case, gn = Writing Cj = Cj^^j^^^i we then have 

fc-i 

1=0 

On multiplying by 2'^^^~^^~^ this yields the equation 

i=0 

with integer entries Cj. We now check that this equation is impossible (mod 2). All terms 
except the last are even since {m — l){k — 1 — i) > when < z < A; — 2. The {k — l)-st term in 
the last sum is Cjfc_i2° = 1, so we get = 1 (mod 2), the desired contradiction. We conclude 
that S 7^ In+2) and the theorem follows. ■ 

6. 5-Integral Apollonian Cluster Ensembles (Dimension n) 

We study integrality and rationality properties for Apollonian cluster ensembles. Given a finite 

set of primes S, we say that a rational number is S -integral if its denominator is divisible only 
by powers of primes in S. By convention we let .S = 1 denote the case when there are no 
primes in S. 

The Apollonian group consists of integer matrices in dimensions 2 and 3, and retains an 
S'-integral structure in all dimensions, for suitable S. In each dimension n we consider the 
questions: 

(1) Does there exist some S and an Apollonian cluster ensemble all of whose Descartes 
configurations consist of spheres having S-integral curvatures ("S'-integral ensemble")? 

(2) Does there exist some S and an Apollonian cluster ensemble all of whose Descartes 
configurations have augmented curvature-center coordinate matrices Wj? S'-integral ("super 
S'-integral ensemble")? 

We show that S-integral ensembles exist in all dimensions, if S is chosen properly. However 
we show that super S-integral ensembles can exist only in dimensions n = 2m^ orn = (2m+l)^, 
for integer m, again with S chosen appropriately. 



19 



6.1. S'-integral Apollonian Cluster Ensembles 



We say that an Apollonian cluster ensemble is S-integral if the curvature of every sphere in 
the ensemble is S'-integral. 

Theorem 6.1. In each dimension n > 2 there exists an S-integral Apollonian cluster ensem- 
ble in which S is specified as: 

(1) S is the set of primes dividing n — 1 if n is even. 

(2) S is the set of primes dividing if n is odd. 

Proof. It suffices to show that the Descartes equation 



has a non-zero 5-integral solution b for each n > 2. There is such a configuration D which 
is not only 5- integral, but integral, with curvatures (0, 0, 1, 1, 1). It consists of two parallel 
hyperplanes separated by distance 2 together with n unit spheres whose centers comprise the 
vertices an (n — l)-dimensional simplex in a hyperplane parallel to the two hyperplanes in the 
configuration, and lying midway between them. 

The other Descartes configurations in the Apollonian cluster ensemble and super- Apollonian 
cluster ensemble generated by this configuration are 5-integral, where S is the set of primes 
dividing the denominator of , since they have associated matrices U Wx) for some U in the 
Apollonian group. ■ 

In dimension 2 we can take S = 1, as we saw in ^5]. In that case an Apollonian cluster 
ensemble consists of the Descartes configurations in an Apollonian circle packing. In [JHI we 
studied in various number-theoretic questions related to the integer curvatures appearing in 
integer Apollonian circle packings. 

6.2. Super S'-integral Apollonian Cluster Ensembles 

We say that a Descartes configuration T> is super S-integral if its augmented curvature- 
center coordinate matrix Wxi is S-integral. Similarly, an Apollonian cluster ensemble is super 
S-integral if every Descartes configuration V in the packing has an S'-integral augmented 
curvature-center matrix. The next lemma reduces the question of super S-integrality of an 
Apollonian cluster ensemble to that of a single Descartes configuration. 

Theorem 6.2. // a single Descartes configuration is super S-integral, then the Apollonian 
cluster ensemble it generates is super S'-integral, where S' consists of S together with all 
primes dividing the denominator of . 

Proof. This follows from the fact that the Apollonian group consists of rational matrices whose 
entries have denominators that are divisible only by primes dividing the denominator of 
when put in lowest terms. ■ 

It seems to be a difficult problem to determine for specific S for which dimensions there 
might exist a super 5-integral Descartes configuration. We consider the weaker question of 



QD,nb = 
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whether in a given dimension there exists a super S'-integral Descartes configurations for some 
S. This is the same as the existence of Descartes configurations T> having a rational ACC- 
matrix Wx), and we caU such Descartes configurations super-rational. 

According to the augmented Euchdean Descartes Theorem 13.11 super-rational Descartes 
configurations occur exactly in those dimensions n in which there exists an invertible rational 
matrix W such that 

(6.2) 

that is, the quadratic form Q£),n is rationally equivalent to the form Qw,n We use this fact to 
determine in which dimensions such configurations exist. 

Theorem 6.3. A necessary and sufficient condition on the dimension n for a super S-integral 
Descartes configuration to exist for some S is thatn = 2P or {2k -if for some positive integer 
k. 

To establish this result, we proceed in a series of lemmas. 

Lemma 6.4. Given a Descartes configuration D in M" its associated augmented matrix Wx> 
has 

det(W5?)2 = n2"+l (6.3) 

Proof. This follows from taking determinants in ()3.12|) . since the right side has determinant 
—2"+'' while the left side has determinant det(Wx>)^ det(Qx>,n) and 

detiQv.n) = --. (6.4) 
n 

To verify this last statement, we apply the following row operations to the matrix Q.„. Add 
rows 2 through n + 2 to the first row, to get a new first row that has all entries — ^ . Then add 
this row multiplied by — ^ to each of the other rows. Aside from the first row, the first column 
is zero, and the lower right (n + 1) x (n + 1) matrix is the identity. But this matrix obviously 
has determinant — -. ■ 

n 

Lemma 6.5. // a super-rational Descartes configuration exists in dimension n, then neces- 
sarily n = 2k^ or (2k — 1)^ for some positive integer k. 

Proof. A necessary condition for the existence of a Descartes configuration D whose augmented 
matrix Wx? has rational entries is that det(Wx)) be rational. This requires that n2"+3 be the 
square of a rational number. By Lemma 16.41 this holds for even n if and only if n is twice a 
square, and for odd n if and only if n is an (odd) square. ■ 

To prove the sufficiency of this condition, we use the theory of equivalence of rational 
quadratic forms, cf. Cassels jTj or Conway jH]. We write Q ~(q Q' to mean that the (rational) 
quadratic form Q is rationally equivalent to Q'. To apply the decision procedure, we first 
diagonalize Qn over the rationals, which we do for all n > 2. 



W^Q^ .W = Q 



-4 

2I„ 
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Lemma 6.6. For each n > 2, the Descartes quadratic form QD,n = In+2 — :^ln+2ln+2 ^^■^ 



n — 1 n — 2 2 

QD,n diag{ , 2, 2, 2, -2). 

n n — 1 6 

Proof. We diagonalize the quadratic form as in Conway [HI pp- 92-94]. Set 

M("+2) := Q^^„ = (xo + yo)In+2 - yoln+2ln+2 , 

where xq = , Uo = ^ - At the j-th stage of reduction we will have 

QD,n diag(di,d2, . . . ,d,,M("+2-^-)) , 

where 

M("+2-i) = (a; . + yj.)In+2-i " l„+2-j 1^+2-^ 

for certain xj, yj. The reduction step is 

(w(^'))^m("+2--'')w(^') = diag(dj+i,M("+i--'')) . 
To specify W^-') we first let Wm(a) be the m x m real matrix 

1 a ■ ■ ■ a 



I 



m— 1 



and we set 



-m+2-j 



(6.5) 



(6.6) 



(6.7) 



(6.8) 



Substituting this in (|6.7|) . its left side yields a matrix with the form of the right side with 

dj+i — Xj , 

and with Xj+i, y^+i given by the recursion 



Vj+i 
Xj+i + Vj+i 



Vj + 



(6.9) 
(6.10) 



Solving this recursion, by induction on j, one obtains 

n-j - I 



Vj 



n-j 
1 



n-j 



< i < n - 2, 
< i < n - 2 . 
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This yields the diagonal elements 



n-j - I 
n-j 



1 < J < n - 3 . 



(6.11) 



with 



We find ^2 = 2:3 = | and 



QD,n diag( 



n- 1 



n 



5 ' ' ' ! o ' 



d2,M 



(4)^ 



M(^) = (X„_2 + yn-2)U - yn-2UA 



For the final step in the reduction we use 

iV^(QA2)A^ = diag(2,2,2,-2), 

with 

N = 



1 


1 


-1 


1 


-1 


1 


-1 




-1 


-1 


1 




-1 


1 


-1 


1 



Qd, 



(6.12) 



1 


-1 


1 


1 


1 


1 




1 


1 


-1 


1 


1 


1 


1 


1 


1 



This completes the reduction. ■ 

Proof of Theorem 16.31 The necessity of n = 2A;^ or {2k — 1)^ was proved in Lemma 16.51 
The sufficiency is equivalent to proving that if n = 2k'^ and n = {2k — 1)"^ then 



QD,n — Q Qw,n ■- 



0-4 
-4 
2I„, 



We begin by noting the rational equivalence 

Qw,n diag(-2, 2, • • • , 2, 2) = diag(-2, 2/„+i) 

via the matrix 



(6.13) 



Wo = - 
" 2 



1 1 
1-10 

2Ir, 



By permuting variables we have Qw,n —Q diag(2, 2, • • • ,2, —2). Thus the theorem is equivalent 
to showing that QD,n is rationally equivalent to diag(2, 2, 2, —2). Lemma 16.61 gives 



n — In — 2 3 

QD,n — Q ( , r, 2, 2, 2, -2) , 

n n — 1 2 
~Q (n(n - 1), (n - l)(n - 2), . . . , 3 • 2, 2, 2, 2, -2) , 

using at the last step a conjugacy by W = diag(n, n — 1, . . . , 2, 1, 1, 1, 1). 



(6.14) 
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The Hasse-Minkowski theorem says that two rational quadratic forms of the same dimen- 
sion are equivalent if and only they have the same signature, the ratio of their determinants is 
a nonzero square, and they are p-adically equivalent for all primes p, cf. Conway p. 96ff]. 
Lemma l6 . 61 shows that the signatures of Q_D,n and diag(2, 2, 2, 2, —2) agree, and the hypoth- 
esis n = 2k'^ or n = {2k — 1)^ is exactly the condition that the ratio of their determinants is a 
square of a rational number, and it remains to check the p-adic invariants. 

The p-adic invariants (Tp(Q) are defined (mod 8), and for a diagonal form 
Q = diag(di, ^2, . . . , c?„), one has 

n 

f^p(Q) = ^(^pidj) (mod 8) . (6.15) 

i=i 

We recall formulas for ap{d) when d G Z, cf. Conway pp. 94-96]. Write d = bp^ with 
{b,p) = 1. For p > 3, and an even power I = 2j, 



'p 

while for an odd power I = 2j + 1 



ap{d) = = 1 (mod 8) , (6.16) 



p (mod 8) if r-] =1, 
ap{d) ^ { Y/ (6.17) 

p + 4 (mod 8) if fM = -1. 

li p = 2 then for an even power / = 2j, 

a2{d) = h (mod 8), (6.18) 

while for an odd power I = 2j + 1, 

( b if 6 = ±1 (mod 8) , 

f72(d) = <^ (6.19) 

[6 + 4 if 6= ±3 (mod 8). 

Now (|6?TH) gives 

n— 3 

(Tp{ClD,n) = ^p((" - •?')(^ - i - 1)) + 3(7^(2) + ap{-2) (mod 8), 

j=0 



while (|6.13|) gives 



n— 3 

crp{Qw,n) = ^ -7^(2) + 3(Tp(2) + fTp(-2) (mod 
j=0 



To show equality of these, it suffices to show that for all p, 

n—3 n—3 



^ap(2) = ^ap((n-j)(n-i-l)) (mod 8) (6.20) 

j=0 j=0 
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holds whenever n = 2k'^ or n = {2k — 1)^. 

Consider first the case that p > 3 is odd. Then each dp (2) = 1, so 

n— 3 

^ap(2) = n-2 (mod 8). (6.21) 
i=o 

Now if p does not divide {n — j){n — j — 1) then CTp{{n — j)[n — j — 1)) = 1. The terms divisible 
by p occur in blocks of two consecutive terms, and we claim that if p divides j then 

<yp{{j + m + (^pUU - 1)) = 2 (mod 8). (6.22) 

Suppose j = bp'', with where {b,p) = 1 and I > 1. If / is even, both terms on the left side of 
1)6. 22(1 are 1 (mod 8) by (|6.16|) . while if / is odd, then if p = 1 (mod 4), then (^^^ = 1; so the 

two terms both have values p (resp. p + 4) according as = 1 (resp. -1), and their sum is 

2p = 2 (mod 8). If p = 3 (mod 4), then (^^^ = so exactly one of (^^^ takes the value 
— 1, and the two terms add up to 2p + 4 = 2 (mod 8). Thus 1)6.22(1 follows. Thus adding up 
the right side of (|6.2Up and grouping terms divisible by p in consecutive pairs gives 

n—3 n—3 

^ap((n-i)(n-j-l)) = ^l = n-2 (mod 8). (6.23) 

3=0 j=Q 

There remains an exceptional case where p\n, in which case n(n — 1) is divisible by p and is 
an un-paired term. Since n = 2k'^ or {2k — 1)^, thus p'\n with / even, hence cjp{n{n — 1)) = 1 
in this case, and (|6.23j) holds. This establishes (|6.2()|) for p > 3. 

Now consider the case p = 2. Certainly (12(2) = 1 so ((6.211) holds. We claim that 

a2{{23 + l)2j) + a2{2j{2j - 1)) = (mod 8) . (6.24) 

Write 2j = 2^6 with b odd, and by checking all possible cases using ((6.18(1 and ((6.19(1 . one 
verifies ((6.24(1 . Suppose n = 2k'^. Then in the right side of ((6.2U(1 all terms pair except the first 
and last, and ((6.24(1 yields 

n—3 

<y2{{n - j{n - j - 1)) = a2{n{n - 1)) +^2(3- 2) 

3=0 

r -1 + -1 if k = (mod 2) 

\ 1 + -1 if /c = 1 (mod 2) 
= n — 2 (mod 8) , 

so ((6.20() holds. If n = {2k — 1)^ = 1 (mod 8) then all term pair except the last term, and 
(HnUl yields 

n—3 

Y,a2{{n-j){n-j-l))=a2{3-2) = -1 (mod 8), 
3=0 
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so 1)6. 20() holds in this case. ■ 

Theorem 16.31 est abhshes the existence of super-rational Descartes configurations in the given 
dimensions, but does not give a bound for the denominators of the rationals appearing in these 
configurations. It could be that in certain dimensions n = 2k'^ and (2A;+1)^ there exist strongly 
integral Descartes configurations, i.e. ones with S = 1. However, even if such configurations 
exist for some n > 2, then the Apollonian cluster ensemble containing them would not inherit 
the super- integrality property, but only super S -integrality as in Lemma l6.2l We leave this as 
an open problem. 

7. n-Dimensional Duality Operator 

In the two-dimensional case we showed that there exists a duality operator D G Aut{Q£)) which 
took each Descartes configuration D to a new Descartes configuration D' that consists of four 
circles orthogonal to the original circles, each passing through three of the intersection points 
of the circles of T). We showed that D was contained in the normalizer of the super- Apollonian 
group in two dimensions. 

The geometric duality operation based on orthogonal spheres generalizes to higher dimen- 
sions as follows. Given n + 1 mutually tangent spheres in n dimensions, there is a unique 
sphere through their points of tangency, and this sphere is orthogonal to each of the given 
n + 1 spheres, as given in the following (known) result. 

Theorem 7.1. Given n + 1 mutually tangent (n — l)-spheres {Cj : 1 < j < n + 1} in M" 

having disjoint interiors, there exists a unique (n — l)-sphere C"*- passing through the "^"^"^"^^ 
tangency points of these spheres. At each such tangency point the normal to the sphere C"*" is 
perpendicular to the normals of the two spheres Ci and Cj tangent there. 

Proof. The assumption of disjoint interiors (we allow interior to be defined as "exterior" 
for one sphere) is equivalent to all ^^^^^^ tangency points of the spheres being distinct. For 
dimension n = 2 there is a unique circle through any three distinct points. However for n > 3 
the conditions are over-determined, since n + 1 distinct points already determine a unique 
(n — l)-sphere, and the main issue is existence. 

Both assertions of the theorem are invariant under Mobius transformations (which preserve 
angles), and there exists a Mobius transformation taking a set of n+1 mutually tangent (n — 1)- 
spheres in having disjoint interiors to any other such set, cf. Wilker j35| Theorem 3]. Thus 
it suffices to prove the result for a single such configuration, and we consider the configuration 
of n + 1 mutually touching spheres of equal radii whose centers are at the vertices of a regular 
n-simplex, and tangency points of the spheres are the midpoints of its edges. The first assertion 
of the theorem holds in this case because there is an (n — 1)- sphere whose center is at the center 
of gravity of this simplex, which passes through the midpoints of every edge of the simplex. 
Indeed the isometrics preserving an n-simplex fix the center of gravity and act transitively on 
the edges. Note that for n = 2 the simplex is an equilateral triangle and C"*- is the inscribed 
circle; however for n > 3 the sphere C"*" is neither inscribed nor circumscribed about this 
simplex. 
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For the second assertion of the theorem, in this configuration the sphere has each edge 
of the n-simplex lying in a tangent plane to the sphere; so the normal to C"*- at the midpoint 
of an edge is perpendicular to that edge. Two spheres C, and Cj intersect at the midpoint of 
an edge, and the normal to their tangent planes points along this edge; thus this normal is 
perpendicular to the normal to C"*- there. ■ 

The second assertion in Theorem 17. II explains why the sphere is termed " orthogonal." 
Thus, given a Descartes configuration of n + 2 spheres Cj, we get a system of n + 2 "orthogonal" 
spheres 

where is associated to the n + 1 spheres obtained by deleting Cj. When n = 2 the new 
spheres are mutually tangent and give a new Descartes configuration; this gives the "duality" 
operation D. For n > 3, however, the spheres are not mutually tangent. In fact for all n 
their curvatures satisfy a relation similar in form to the original (two-dimensional) Descartes 
relation, namely 

n+2 ^ ^ n+2 



r- 2 ^-^ r, 

i=i « i=i • 



and not the Soddy-Gossett relation 

n+2 n+2 

1=1 ' t=i 

satisfied by Descartes configurations in n-dimensions, cf. |231 Theorem 1.2]. (We omit a 
proof of (|7.1() .) In particular, for n > 3 given a Descartes configuration V, the set V-^ := 
{C^, . . . , C^_^_2} of orthogonal spheres is not a Descartes configuration. 

The question arises, are these n + 2 "orthogonal" spheres in any special relation to one 
another? We answer this in terms of an inversive invariant of two arbitrary (not necessarily 
tangent) oriented spheres. 

Definition 7.1. (i) The separation between two oriented spheres Ci and C2 with finite radii 
ri and r2, and with centers distance d apart, is 

A(Ci,C.):=^^^i^, (7.3) 
2rir2 

provided both spheres are inwardly oriented or outwardly oriented, and is otherwise the nega- 
tive of the right side of this formula. 

(ii) The separation of an oriented sphere Ci of finite radius ri and an oriented hyperplane 
C2 is 

A(Ci,C2):=-, (7.4) 

ri 

where d is the (signed) distance from the center xi of Ci to C2, measured so that d > if xi 
is not in the interior of C2 and Ci is inwardly oriented, or if xi is in the interior of C2 and Ci 
is outwardly oriented, and d < otherwise. 
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(iii) The separation between two oriented hyperplanes Ci and C2 is 



A(C7i,C2) := -COS0, (7.5) 

where 9 is the dihedral angle between the designated normals at a point of intersection. 

The separation of two spheres is an inversive invariant, hence an Mobius invariant; that is, 

A(0(C7i),0(C72)) = A(Ci,C2), (7.6) 

holds for any Mobius transformation g. This concept appears in Boyd who introduced 
the term separation for it, but the concept ^ was used earlier by Mauldon j2S] in 1962, who 
used the term inclination to mean the negative of A(Ci,C2), and showed it was an inversive 
invariant. The absolute value of A(Ci,C2) was studied in Beardon ^ pp.29] under the name 
of inversive product of two spheres. 

The separation A(Ci,C2) of two spheres can be expressed in terms of their augmented 
curvature-center coordinates as 

A(Ci,C2) = ^w(C7i)^K„w(C2) 

1 - - " 

MCi)b{C2) + b{Ci)biC2))+b{Ci)biC2)^XjiCi)xjiC2), (7.7) 



2 



where K„(l) is given by 



0-10 
-10 
2/„ 



(7.8) 



This formula can be proved by a simple algebraic calculation. Using it, one can check that 
for two tangent spheres Ci and C2, A(Ci,C2) = 1, if (1) Ci and C2 are externally tangent, 
and both are inwardly oriented or outwardly oriented, or (2) Ci and C2 are internally tangent 
and one is inwardly oriented, the other is outwardly oriented. In all other cases two tangent 
spheres have A(Ci,C2) = —1, and orthogonal spheres are those with A(Ci,C2) = 0. 
Prom Theorem 17. II one obtains 

A(C^, Cj) = for 1 < J < n + 1, (7.9) 

and these relations determine C"*" up to orientation. It can also be shown that if a set of 
tangent spheres {Ci, C^+i} have oriented curvatures b^+i = (61, 6„+i), and centers Xj, 
then for either orientation the orthogonal sphere C"*" has oriented curvature q satisfying 

(7.10) 




■^The idea of considering such an inversive invariant traces back to work of Clifford jgj in 1868 and of Darboux 
|14| in 1872. However, neither Clifford's nor Darboux' definition was precisely A(Ci,C2). Clifi'ord defines the 
power of two spheres to be the square distance of their centers less the sum of the squares of their radii, i.e., 

— — rl, and Darboux also uses the same quantity, ^14! p. 350]. 
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and (oriented) center x satisfying 

gx = -b„+i(^QD,n-i)C, (7.11) 

in which C is an (n + 1) x n matrix whose j-th row is bjXj. 

An oriented Descartes configuration in is characterized in terms of separation as a set 
of n + 2 oriented spheres each pair of which has A{Ci,Cj) = 1, when i ^ j. Thus such a 
configuration has the following property. 

Definition 7.2. A collection of oriented spheres is equiseparated if all values A(Cj, Ck) with 
j ^ k are equal. 

The equiseparation property can also be viewed as an equiangularity property, because for 
two oriented circles that intersect or touch one has 

A(C7i,C72) = -COS0, (7.12) 

where is the angle between oriented normals at a point of intersection of the two circles. We 
now show the duality operation preserves equisepar ability in all dimensions. 

Tlieorem 7.2 (Equiseparation Tlieorem) Given an oriented Descartes configuration D = 
(Ci, C2, C„4.2) in M", if the dual spheres are properly oriented then the (oriented) dual 
configuration (Cf*", C^, C^_^2) ^-^ equiseparated, with 

A{Cf,Ct) = -^ if j^k. (7.13) 

■' n — 1 

Proof. In this result, the orientation assigned to the dual spheres in the theorem depends on all 
n + 2 spheres in the Descartes configuration, and the orientation of Cj' cannot be consistently 
assigned from the n + 1 oriented spheres {Cj ■ i ^ j} alone. If all n + 2 spheres Cj are inwardly 
oriented, then n + 1 of the spheres Cj~ will be inwardly oriented and one outwardly oriented, 
the last being the one of largest radius. If all but one of the n + 2 spheres are inwardly oriented, 
and one outwardly oriented, then all n + 2 spheres Cj*- will be inwardly oriented. 

Since the result is invariant under inversion, it suffices to prove it for a single Descartes 
configuration. We consider the special oriented Descartes configuration where the curvatures 
are (0,0,1,1,. . . , 1). Here we have two parallel planes, which we take as xi = ±1, and n unit 
spheres, all with centers on the plane xi = 0. Their centers form a regular simplex in this 
plane. We may take one of these centers at (0,^,0,0, . . .) where = 2(n — l)/n. Consider 
the "orthogonal" spheres that pass through the point T = (1, ^, 0, 0, . . . , 0). There are n such 
spheres, and all but one of them is a hyperplane containing the points T and (—1, ^, 0, . . . , 0), 
and the centers of all but one of the original unit spheres. These centers form the vertices of a 
regular (n — l)-simplex, so these n — 1 "orthogonal" planes are equiangular satisfying H7.12() . 
where 9 is the angle between the normals of two facets of a regular n-simplex. It follows that 
these orthogonal planes satisfy (|7.13|) . The final "orthogonal" sphere through T is orthogonal 
to the plane xi = 1 and all the n original unit spheres. Its center is thus (1,0,0, . . .) and its 
radius is ^. Hence it is also equiangular with the n — 1 "orthogonal" planes, with cos 9n = 
These angles are all equal to the one formed by connecting the vertices of a regular simplex 
to its center, i.e. the angle in a triangle of sides ^, ^ and 2. Finally, the last two "orthogonal" 
spheres meet at the same angle in the plane xi = 0. ■ 
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8. Concluding remarks 



This paper studied generalizations of the basic properties of two-dimensional Apollonian pack- 
ings and super-Apollonian packings. We gave fairly complete answers but left open a few 
problems. 

The first problem is to determine a presentation of the super-Apollonian group in 
dimensions n > 3, in terms of the given generators. Is this group always a hyperbolic Coxeter 
group? 

There are some unanswered questions concerning rational and integral structures on Descartes 
configurations. In §3.1 we raised the question of determining those dimensions n in which the 
Descartes form QD,n, the Wilker form Qw,n and the Lorcntz form Qc,n are rationally equiv- 
alent. In §6 we gave a necessary and sufficient condition for rational equivalence of the pair 
(Descartes, Wilker), namely that n = 2k^ or n = (2fc — 1)^. For the other two pairs (Descartes, 
Lorentz) and (Wilker, Lorentz), we gave the necessary conditions n = 2k'^ and n = 2k, respec- 
tively. It remains to determine necessary and sufficient conditions in these cases. 

Finally we left open the question of whether there is any dimension n > 3 in which there 
exist strongly integral Descartes configurations. Although strong integrality will not be pre- 
served under the action of the super- Apollonian group (since n > 4), it would be preserved 
under the action of the dual Apollonian group A^, which consists of integer matrices in all 
dimensions. 
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9. Appendix. Mobius Group Action 



The Mobius group action given in Appendix A of part I straightforwardly extends to n- 
dimensions. 

The (general) Mobius group Mob(n) is the group generated by reflections in spheres or 
planes in the one-point compactification = U {00} of M". see Beardon (21 Chapter 
3]. (He denotes it GM(M„).) This group has two connected components, and we let Mob(n)_(_ 
denote the connected component of the identity. The extended Mobius group GM*{n) is defined 
by GM*{n) := Mob(n) x {—1,1}. Here {—1,1} are in the center of this group, and we write 
elements of GM*{n) as itg, in which g G Mob(n), and the sign indicates which of ±1 occurs. 
The group GM*{n) has four connected components. 

The purpose of this Appendix is to define an action of GM* (n) on the right on the parameter 
space A^jj, given in Theorem 19.11 below. This amounts to finding an explicit isomorphism 
between GM*{n) and Aut{Qw,n)- The case n = 2 was treated in Appendix A of part I ( |18| 
Theorem 7.2]), was shown that the Mob(2) action preserves (total) orientation of Descartes 
configurations. The same property holds for the Mob(n) action treated here, by a similar proof 
which we omit. 

Relevant isomorphisms are pictured in Tabled below. The isomorphism between GM*{n) 
and Aut{Qw,n) appears as the horizontal arrow on the left in the top row. This map when 
restricted to the smaller groups Mob(n) and Mob(n)+ give the other two horizontal isomor- 
phisms on the left side of the table. Table H also indicates isomorphisms on its right side to 
the orthogonal group 0(n + 1, 1) and corresponding subgroups, which we defer discussing until 
after the following result. 

Theorem 9.1. Let GM*{n) := Mo6(n)x{I, —I}. There is a unique isomorphism vr : GM*{n) - 
Aut{Qw,n) , with image elements Vj-g := 7r(ibg), such that the following hold. 

(i) For Q S Mdb{n) the augmented curvature- center coordinates for each ordered, oriented 
Descartes configuration V satisfy 

W0(^) =Wi,Vgi. (9.1) 

(a) The action of —I on augmented curvature- center coordinates is 

W_75 = WpVl} = -Wt). (9.2) 
This reverses the orientation of the Descartes configuration. 

Proof. We compute the action of Mob(n) on augmented curvature-center coordinates. Let 

{b, b, wi,W2, . . . , Wn) = ( , ^7) 77, • • • , ^) be the augmented curvature-center co- 
ordinates of the sphere 

n 

^{Vi - Xif = 

i=l 

This sphere can be recovered from these coordinates via 

n 

Y^ihyi - Wi)^ = 1, (9.3) 

2 = 1 
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and the orientation of the sphere (inside versus outside) is determined by the sign of b. An 
oriented "sphere at infinity" is a hyperplane given by 



y • h = m, 

and its associated curvature-center coordinates are 

{b, b,Wi,W2,...,Wn) = (2171,0, hi, h2,...,hn), 



(9.4) 



(9.5) 



where h = {hi,h2, ■ ■ ■ ,hn) is the unit normal vector, and the orientation is given by the 
convention that the normal h points inward. 
The group M6b(n) is generated by 

(1) translations tyo(y) = y + Yo; 

(2) dilations X)r{y) = ry with r G M, r > 0; 

(3) the rotation o(y) = Oy, where O is an n x n orthogonal matrix; 

(4) the inversion in the unit circle jc(y) = j^- 

Given g G M6b(n), we will let § denote the corresponding action on the curvature-center 
coordinates of an oriented circle. The action of translation by yo = (2/0,1)2/0,2, • • • ,yo,n) is 



tyg{b, b, Wi,W2,...,Wn) 



{b + 2'^Wiyo^i + b'^yQ^i,b,wi + byo,i,W2 + 61/0,2, • • • , w^n + 6yo,n)- 



i=l 



i=l 



The action of a dilation with r G M, (r > 0) is given by 

X)x(b,b,wi,W2, . . .,Wn) = {rb,b/r, wi,W2, . . .,Wn)- 
The action of rotation with orthogonal matrix O is 

o{b,b,wi,W2,...,Wn) = {b,b,w[,w'2,...,w'J, 
where {w[,W2, . . . ,w'^) = {wi,W2, ■ ■ ■ , Wn)0'^. The action of inversion in the unit circle is 

jcib,b,Wi,W2, . . .,Wn) = {b,b,'Wi,W2, . . .,Wn). 

All of these actions apply to "spheres at infinity" and extend to linear maps on the (n + 2) x 
(n + 2) matrices Wx>. 

The translation operation is given by right multiplication by the matrix 



V 



-1 



1 

i=iyo,i 1 2/0,1 

22/0,1 •■• 

2yo,n 
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and one verifies H9.1|) holds by direct computation. For the dilation i)^, with r G M, r > the 
right action is by the matrix 



r 
1/r 
I„ 



For rotation o, the right action is by the matrix 



* 



1 
1 




For the inversion jc the unit circle, the permutation matrix 



V 



V 



(12) 



1 

1 
I„ 



It is easy to verify that the above matrices are all in Aut{Q\Y). Since Q,w,n 
where 

"110 
A = V2 -11 
. I„ 

we have that a matrix U G Aut{Qw,n)^ if and only if e^AUA^-'^ei > 0. One checks that all 
the above matrices are actually in Aut{Qw,nV ^ so that the map so far defines a homomorphism 
of Mob(n) into Aut{Qw^n)^ — 0{n + l, 1^ , identified with the isochronous Lorentz group. The 
group Mob(n) acts simply transitively on ordered Descartes configurations, as observed by 
Wilker [HEl Theorem 3, p. 394], and the group Aut{Q]Y^n) acts simply transitively on ordered, 
oriented Descartes configurations, as implied by Theorem 3.1. Because Aut{Qw,n)^ is of index 
2 in Aut{Qw,n) — 0{n + 1, 1), we conclude that the map so far defines an isomorphism of 
Mob(n) onto Aut{Q 



W,n) 



To complete the proof, we define the action of — / to be 



,-1 



*-n+2 ■ 



(9.6) 



It has the effect of reversing (total) orientation of the Descartes configuration, and does not 
correspond to a conformal transformation. Since — In+2 ^ Aut{Qw,n)\ adding it gives the 
desired isomorphism of GM*{n) onto Aut{Q\\rn). ■ 



We now return to the data in Table giving the isomorphisms of Aut{Qw,n) and its 
subgroups to the orthogonal group 0{n + 1,1) and its two subgroups 0(n + 1,1)^ the or- 
thochronous orthogonal group, and 0{n + 1, 1)^ the special orthochronous orthogonal group, 
which is the connected component of the identity of the orthogonal group 0(n + 1, 1). The 
set of isomorphisms given by the three horizontal arrows on the right in Tableware obtained 
by any fixed choice of real matrix A that intertwines Qw,n and Qc,n by Qvy = A-^Q£_„A, in 
which case the isomorphism is Aut{Qy/^n) = A~^0{n + 1,1) A sending V i— > AVA~^. Such 
matrices A exist in all dimensions. It is shown in §3.1 that rational matrices A do not exist 
in all dimensions; a necessary condition for existence is that the dimension n be even. 
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GM*(n) = Mdb{n) x {-1,1} — - AutiQw,n) 0{n + 1, 1) 



TT 



TT 



TT 



Mdh{n) 



Aut{Qw,n)^ ^ 0{n + 1, 1) 



Mdb{n)^ 



Aut{Qw,n)l ^ 0{n + 1, 1)1 



Orthogonal group 



Orthochronous 
orthogonal group 



Special orthochronous 
orthogonal group 



Table 1: Group Isomorphisms 



34 



References 

[1] D. Aharonov and K. Stephenson, Geometric sequences of discs in the Apollonian packing. 
Algebra i Analiz 9 (1997), No. 3, 104-140. [Enghsh version: St. Petersburg Math. J. 9 
(1998), 509-545.] 

[2] A. F. Beardon, The Geometry of Discrete Groups, Springer- Verlag: New York 1983. 
[3] M. Berger, Geometry II, Springer- Verlag: Berlin 1987. 

[4] D. W. Boyd, The osculatory packing of a three-dimensional sphere. Canadian J. Math. 
25 (1973), 303-322. 

[5] D. W. Boyd, A new class of infinite sphere packings. Pacific J. Math. 50 (1974), 383-398. 

[6] R. W. Brooks, The spectral geometry of the Apollonian packing. Comm. Pure. Appl. 
Math. 38 (1985), 359-366. 

[7] J. W. S. Cassels, Rational Quadratic Forms, Academic Press: New York 1978. 

[8] W. K. Clifford, On the powers of spheres (1868), in: Mathematical Papers of William 
Kingdon Clifford, MacMillan and Co., London 1882, pp. 332-336. 

[9] J. H. Conway, with F. Fung, The sensual (quadratic) form, Carus Monograph No. 26, 
Math. Assoc. America, Washington DC, 1997. 

[10] H. S. M. Coxeter, The problem of Apollonius. Amer. Math. Monthly 75 (1968), 5-15. 

[11] H. S. M. Coxeter, Introduction to Geometry, Second Edition, John Wiley and Sons, New 
York, 1969. 

[12] H. S. M. Coxeter, Loxodromic sequences of tangent spheres. Aequationes Mathematicae 
1 (1968), 104-121. 

[13] H. S. M. Coxeter, Numerical distances among the spheres in a loxodromic sequence. The 
Mathematical Intelligencer 19 (1997), 41-47. 

[14] G. Darboux, Sur les relations entre les groupes de points, de cercles et de spheres dans le 
plan et dans I'espace, Ann. Sci. Ecole Norm. Sup. 1 (1872), 323-392. 

[15] K. J. Falconer, The Geometry of Fractal Sets. Cambridge Tracts in Math., vol. 85, Camb. 
Univ. Press, Cambridge, 1986. 

[16] T. Gossett, The Kiss Precise, Nature 139(1937), 62. 

[17] T. Gossett, The Hexlet, Nature 139 (1937), 251. 

[18] R. L. Graham, J. C. Lagarias, C. L. Mallows, A. Wilks and C. Yan, Apollonian circle pack- 
ings: geometry and group theory I. Apollonian group, eprint: arXiv lmatIi.MG/ 00 102981 



35 



[19] R. L. Graham, J. C. Lagarias, C. L. Mallows, A. Wilks and C. Yan, Apollonian circle 
packings: geometry and group theory II. Super-Apollonian group and integral packings, 
eprint: arXiv |niath.MG/0010302| 

[20] R. L. Graham, J. C. Lagarias, C. L. Mallows, A. Wilks and C. Yan, Apollonian cir- 
cle packings: number theory, J. Number Theory 100 (2003), 1-45. [eprint: arXiv 
niath.MT/009113] 

[21] E. Kasner and F. Supnick, The Apollonian packing of circles. Proc. Nat. Acad. Sci. USA 
29 (1943), 378-384. 

[22] R. Lachlan, On systems of circles and spheres, Phil. Trans. Roy. Soc. London, Ser. A 177 
(1886), 481-625. 

[23] J. C. Lagarias, C. L. Mallows and A. Wilks, Beyond the Descartes circle theorem, Amer. 
Math. Monthly 109 (2002), 338-361. [eprint: arXiv math. M G/0 10 10661 

[24] B. B. Mandelbrot, The Fractal Geometry of Nature, Freeman: New York, 1982. 

[25] J.G. Mauldon, Sets of equally inclined spheres, Canadian J. Math. 14 (1962), 509-516. 

[26] G. Maxwell, Sphere packings and hyperbolic reflection groups. J. Algebra 79 (1982), 78- 
97. 

[27] Z. A. Melzak, Infinite packings of disks. Canad. J. Math. 18 (1966), 838-853. 

[28] D. Pedoe, On a theorem in geometry, Amer. Math. Monthly 74 (1967), 627-640. 

[29] B. Rodin and D. Sullivan, The convergence of circle packings to the Riemann mapping. 
J. Differential Geometry 26 (1987), 349-360. 

[30] F. Soddy, The Kiss Precise. Nature 137 (1936), 1021. 

[31] F. Soddy, The Hexlet, Nature 138 (1936), 958. 

[32] F. Soddy, The Bowl of Integers and the Hexlet, Nature 139 (1937), 77-79. 
[33] F. Soddy, The Hexlet, Nature 139 (1937), 154. 

[34] B. Soderberg, Apollonian tiling, the Lorentz group, and regular trees. Phys. Rev. A 46 
(1992), No. 4, 1859-1866. 

[35] J. B. Wilker, Inversive Geometry, in: The Geometric Vein, (C. Davis, B. Griinbaum, F. 
A. Sherk, Eds.), Springer- Verlag: New York 1981, pp. 379-442. 

email: graham@ucsd.edu 

j agariasQumich . edu 
colinmSresearch . avayalabs . com 
allanSresearch . att . com 
cyanSmath . tamu . edu 



36 



